function [act logact] = gmmactivRobust(mix, X)


K = mix.ncentres;
N = size(X, 1);

act = zeros(N, K);
logact = zeros(N, K);
for k = 1:mix.ncentres
    [act(:, k) logact(:, k)] = gaussRobust(mix.centres(k, :), mix.covars(:, :, k), X);
end
